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Abstract. Polymer adsorption on fractally rough walls of varying dimensionality is 
studied by renormalization group methods on hierarchical lattices. Exact results are 
obtained for deterministic walls. The adsorption transition is found continuous for 
low dimension d w of the adsorbing wall and the corresponding crossover exponent 
rf> monotonically increases with d w , eventually overcoming previously conjectured 
bounds. For d w exceeding a threshold value d^, <j> becomes 1 and the transition 
turns first-order, d^ > d saw , the fractal dimension of the polymer in the bulk. An 
accurate numerical approach to the same problem with random walls gives evidence 
of the same scenario. 
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1. Introduction 



The adsorption on an attracting impenetrable wall is perhaps the most elementary 
transition involving a single interacting polymer in solution Q. High dilution in 
a good solvent is the realistic condition for which this problem can be directly 
relevant. The fundamental character and the obvious relation with more complex 
applications, like colloid stabilization or surface protection |^|, attracted on polymer 
adsorption a great deal of attention in recent years, and much information is presently 
available on this problem. It is now well understood that this transition can be 
interpreted as a surface critical phenomenon jl], at the adsorption temperature 
T a the conformational statistics of the polymer shows a multicritical behavior with 
peculiar geometric features and with crossovers to the high-T desorbed and the low- 
T adsorbed regimes. For a chain with N monomers at T a the average number of 
adsorbed monomers, (M), scales as (M) oc N^, where <fi (0 < cj) < 1) is the crossover 
exponent. In the high- and low-T regimes, (M) oc iV° and (M) cx N, respectively. 
4> is known exactly in 2D for a polymer in both good and theta || solvents, and 
in 3D in theta solvent, in which case logarithmic corrections are present |J. Further 
exact results have been obtained for models defined on fractal lattices, like Sierpinski 
gaskets [0, ^, § , which are by now recognized as an important context in which to test 
theoretical ideas concerning polymer statistics. 

Most explicit results obtained so far on polymer adsorption refer to cases in which 
the wall is smooth and flat. In this paper we address the adsorption transition on a 
fractal substrate. This problem has applicative interest. Indeed, in many processes 
involving polymers, highly corrugated, irregular walls may be present. In addition 
there are interesting theoretical implications. A polymer in good solvent is known to 
possess a self-similar stochastic geometry, with a well defined fractal dimension d saw ■ 
Once such a polymer is put in contact with a fractal wall, a competition between 
the two geometries arises. This holds in particular for cases in 2D or in hierarchical 
lattices, for which the wall is topologically one dimensional like the polymer. As we 
will see, the above competition leads to a modification of the universal properties of the 
adsorption transition, or, in more extreme situations, to a drastic and quite unusual 
suppression of its continuous character. The parameter triggering such modifications 
is the fractal dimension d w of the wall. This scenario has also analogies with another 
situation in which competition between two similar scaling geometries has been studied 
recently. A fluctuating interface between two coexisting phases is self-affinc Jh]] . If it 
is put in contact with a rough wall of similar geometry, depinning from the wall turns 
from continuous to first-order as soon as the roughness exponent of the wall exceeds 
the anisotropy index of interface fluctuations in the bulk jll], 

An approach to our problem on Euclidean lattice would meet very serious 
difficulties. Once assigned a given profile to a fractal wall exerting short range 
attraction on a self-avoiding chain (SAW), exact enumerations would handle too 
short chains, unable to feel the fractal corrugations of the wall on a sufficiently wide 
range of length scales. On the other hand, Monte Carlo simulations meet the serious 
obstacle that, in the low-T region, sampling over polymer configurations becomes very 
problematic, due to the highly irregular wall, with valleys and hills at all scales. For 
these reasons adsorption on fractal walls is certainly one of those phenomena for which 
the study on simplified, hierarchical lattices is at present the only realistic way to gain 
an at least qualitative understanding. Two recent works studied adsorption on fractal 



boundaries of SAW's within fractal lattices 13, 14 . Most emphasis there was put 
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on the existence of violations of bounds suggested in previous work [j7| for (f>. It was 
also realized that fractal lattices with peculiar connectivities at the borders could give 
rise to an interesting dependence of <fr on the interaction parameters. However, such 
nonuniversality, met also with flat walls, is specific of the lattices considered, which 
do not mimic generic situations. Indeed, it is natural to expect universal scaling at 
the adsorption transition for a given (universal) bulk criticality and a given boundary 
geometry (no matter whether flat or rough). Changes of <p should be expected upon 
varying d w . This is indeed the dimension pertaining to the surface critical phenomenon 
to which adsorption amounts. 

Here we study adsorption in three hierarchical lattices leading to renormalization 
group (RG) recursions of increasing complexity and, supposedly, to results of 
increasing qualitative value with reference to realistic situations on 2D Euclidean 
lattice. 



2. Models and Results 

Let us consider first the lattice Ca, whose construction rule is sketched in Figure |l|a. 
Measuring the lattice "linear size" in terms of the number of steps of the shortest path 
between top and bottom vertices, at any application of the construction rule, this size 
and the total number of lattice bonds are multiplied by factors 2 and 5, respectively. 
Thus, Ca has a fractal dimension dh = In 5/ In 2 = 2.322 ... At any level, n, in the 
construction of Ca, allowed polymer configurations correspond to SAW's between the 
top and bottom vertices. In a grand canonical formulation, to each step is associated 
a monomer fugacity u = exp(///T). An attracting impenetrable wall is modeled as a 
particular SAW which can not be trespassed by the polymer. The polymer interacts 
with the wall through an attractive contact potential — e. Thus, SAW steps on the 
wall acquire an extra fugacity k — exp(e/T) (k > 1). 

In the absence of wall, through the n — lattice there is a unique walk of unit 
length and the restricted grand partition function for SAW's joining top to bottom 
is simply Bq = oj. At level n — 1 there are two pairs of SAW's of lengths 2 and 3, 
respectively. The corresponding partition function is B\ — 2uj 2 + 2ui 3 = 2B 2 + 2Bq = 
B(Bq). If we denote by B n the SAW partition at level n, we can write B n+ i = B(B n ) 
and B can be seen as a generating function of the bulk partition function. The 
recursion for B n has a repulsive fixed point at B* = (V3 — l)/2 = 0.366 . . . which 
corresponds to the bulk critical point of SAW's. Thus u c = B* is the SAW critical 
fugacity. For an n-th level lattice the average number of SAW steps is given by 
(N B )n = jr lit- At w = w c we have (N B ) n = X% with X B = ^jp- = 2.268 . . . 
Therefore, taking into account that the lattice size is L n — 2™, we conclude that 
critical SAW's are fractal with dimension d saw = lnAs/ln2 = 1.181 . . . 

Now let us consider a wall through Ca- In Figure lb and lc we sketch two 
examples of deterministic rules by which wall geometries with opposite features can 
be realized. Iteration of the rule sketched in Figure lb produces a wall whose length 
is L n . Thus, this wall is characterized by a dimension d w = 1 and we regard it as 
fiat. On the contrary, the rule in Figure lc produces a fractal wall with dimension 
d w = In3/ln2 = 1.585... which is also the highest realizable in Ca- Walls with 
intermediate dimensions can be obtained by using either deterministic, or random 
sequences of the two rules above at progressing levels of lattice construction. As 
an example, let us consider a case in which the wall is realized by means of rule 
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lc and lb for odd and even n, respectively. The resulting wall has a dimension 
d w = In 6/ In 4 = 1.292 . . . and is sketched in Figure Id. For n = 1 the SAW partition 
function in the presence of the wall is given by 



Xy — Xq + XqBq — X2(Bq 7 Xq) 

where Xq — ku; is the partition of the unique SAW in the n 
have 



lattice. At 



X 2 = Xf + 2XyB{ 



By = Xy{By,Xy). 



(1) 

2 we 



(2) 



Xy and X2 are now the generating functions of the SAW partitions corresponding to 
rules lb and lc, respectively. At any other construction level, the form of the recursions 
is the same as in Equation (Q) or Equation (0) for n even or odd, respectively. Focusing 
attention on even n, if X n denotes the SAW partition at level n, we have 



X,, 



X 2 (B(B n ),Xy(B n ,X n )). 



(3) 



For B = uj c the bulk recursion is at its repulsive fixed point, while (||) has an 
attractive fixed point at X = 0.1536 . . . and a repulsive one at X = x c = 0.7249 . . . 

Another attractive fixed point is at X = 00. While the first fixed point controls 
the SAW ordinary desorbed regime, the second one corresponds to the adsorption 
critical point which is then located at a wall attraction k c = x c /oj c = 1.9806 . . . The 
critical exponents can be obtained by linearization of the RG flow around the fixed 
point (uj c ,x c ). In doing this, together with the recursion (^) we have to consider also 
B n+ 2 — B(B(B n )) and the matrix 
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and a — 
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OB„ 



. If we put (N) n = and (M)„ = jr^r, at the transition, 

the average number of steps on the wall grows as (M) n = A". Thus, (M) n ~ L^ 1 with 
yi = In A x / In 2 = 0.8001 . . . On the basis of (0) and (0) one can also write: 



\n — 2 \n — 2 

(N) n = \? + a\l \ 2 Al 2 

Ag - Ay 



(6) 



This relation shows that, as long as Ai < As, in the limit of large n, (N) n scales 
as A^. In other terms (N) n ~ with y = ys- On the other hand, if Ai > A^, 
{N) n ~ Ay and y = lnAi/ln2 = yy. In this latter case, which, as we will see below, 
is never realized for La, one must find yy = d w because (M) n ~ (N) n , which again 
means (f> = 1. At the adsorption fixed point of (||) (k — k c ), Xy < Xb- thus, there, the 
SAW length in the presence of wall scales as in the bulk without wall. The previous 
results allow to write (M) n ~ (N)fc with <j> = yyjys = 0.6773 . . . 

For any k > k c the SAW's are adsorbed on the wall and the scalings of both (M) n 
and {N) n are controlled by a fixed point at infinity. In this case both quantities are 
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proportional to . On the other hand, for k < k c , we are in the normal regime 
in which (N) n ~ L^ B . For SAW's on Euclidean lattice, in this last regime, (M) n 
saturates to a constant value. To the contrary, the analysis of recursion (|^) around 
X = 0.1536 . . . shows that (M) n ~ L~ 1 0844 ". This unphysical behavior is due to the 
pathological increase of coordination with n — > oo typical of hierarchical lattices. The 
above discussion can of course be adapted to the case of odd n levels, with the same 
scaling results. 

The method just illustrated can be easily generalized to fractal walls with a whole 
range of d w . In fact, adsorption on walls with different dimensions can be obtained 
by simply changing the form of the recursion. In general, in place of (J3[) , we can have 
p and q nested applications of the generating functions X\ and X 2 , respectively. For 
example, for p — 2 and g = 3we could have 

X-n+p+q = 

Xt{B{B{B{B{B n )))\ X 2 (B(B(B(B n ))), X 2 (B(B(B n )), X x {B{B n ), X 2 (B n , *„))))) (7) 

which corresponds to applying successively X 2l X\, X 2 , X 2l and X±. The corresponding 
wall has d w — (^^ni^j ■ F° r given p and q, we could generate many different 
recursions by changing the order in which X\ and X 2 are applied. While this order 
influences the position of the adsorption transition point, the corresponding critical 
exponents are only functions of the wall dimension. This property follows from the 
circumstance that by construction the different RG recursions are diffeomorphically 
related. This is the mechanism by which d w alone, and no other features of the wall, 
determines the critical exponents of the transition. In these more general conditions 
we find 
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where x c has to be calculated for the specific form of X n+p+q . As already mentioned, 
d w can be varied between 1 and In 3/ In 2 in Ca- The lower limit corresponds to p = 1 
and q = 0. In this case we have k c = 1 and Ai = 1 which implies <f> — 0. In fact, 
for this case of "flat" boundary, the adsorption transition fixed point merges with 
the desorbed regime fixed point and becomes marginally unstable. This is a peculiar 
feature due to the relatively too simple structure of Ca- On the other hand, the upper 
limit corresponds to p = and q = 1 for which one has k c = y 3 + \/3 and Ai = Xb- 
In this case we have exactly 0=1. For intermediate wall dimensions, <f> is a monotonic 
increasing function of d w . A summary of our results for Ca is reported in Table |l|. 

4> < 1 implies continuity of the adsorption transition. On the contrary, for 0=1, 
at the transition (M) n oc (N) n , as for an adsorbed polymer. In fact Ai > Xb implies 
a discontinuity of lim„^ 00 [(M) n /(iV) rl ] at the adsorption point, i.e. a first-order 
transition. This discontinuity, found only for d w = In 3/ In 2 in Ca, anticipates a 
more general result, valid for the other lattices we considered: we find below that a 
sufficiently high d w can drive polymer adsorption first-order. For Ca the threshold 
condition for discontinuous adsorption occurs precisely when d w is at its maximum 
possible value. In general we will denote by d^ the value of d w above which 0=1. 

To show that the scenario described above is not just peculiar of Ca and to 
investigate further the change of transition order, we considered SAW adsorption on 
fractal walls also with Cb and Cc, sketched in Figures ||a and ||a, respectively. Cb has 
a diamond structure similar to that of Ca, but with a higher ramification and d^ = 3. 



() 



The bulk SAW partition function obeys the recursion B n+ i = 3B 2 + 4B 3 + 2B 4 
with lu c = B* — 0.2441 . . . and d saw — 1.1995 ... In Figure ^Jd and |^c we report two 
construction rules generating walls with dimensions dffl — 1 and dffl — 2, respectively. 
The corresponding generating functions are: 

Xi (X, B)=X 2 + 2XB 2 + B 2 (9) 
X 2 (X,B) = XB + X 2 B + X 4 . (10) 

For fractal walls obtained by suitably alternating the two rules, Table || shows again 
that <j) monotonically increases with d w . However, now for any d w > 1.67 we have 
= 1. Thus, when d w > d* w ~ 1.67 adsorption becomes discontinuous. 

Finally, we consider C c for which d L = In 12/ In 4 = 1.792... [Q. The bulk 
recursion is now B n+ i = 6B 4 + 4B® + 2B\ with critical fixed point at lo c = B* = 
0.5175... and d saw — 1.0649... Unfortunately, neither w c , nor d saw are too close 
to the values appropriate for 2D square lattice, 0.378 . . . and 4/3 |l7], ||, respectively. 
This occurs in spite of the fact that Cc seems to mimic well the square lattice structure 
at local level. Cc offers more possibilities of wall construction rules. In Figures ||b, ||c 
and |^d we show three examples with respective generating functions 

Xt(X,B) = 2B 4 + X 4 + 2X 2 B 2 + 2X 2 B 4 (11) 
X 2 (X, B) = XB 3 + X 2 B 4 + X b B 3 + X 2 B 2 + X 3 B + X 6 (12) 
X 3 (X,B) = X 3 B + X 5 B + X 8 (13) 

The corresponding walls have dw^ — 1, d!ff — In6/ln4 = 1.292... and dffl — 
In 8/ In 4 = 1.5. Results for Cc are reported in Table [| For d w = 1 adsorption is 
continuous and <\> — 0.5437 ... is not too far from 1/2, the <j> value for SAW adsorption 
on a smooth wall in 2D Ji] . For increasing d w , cj> monotonically increases and reaches 
a unit value at d w = ~ 1.4. For d w > d^, 4> = 1 an( i the transition is always 
first-order. Note that the first-order transitions found for d w > d^ in Tables 2 and 3 
correspond to Ai > As. From Tables 0, || and || we also learn that, upon increasing 
d w , k c increases, as a rule, up to small fluctuations caused by our peculiar recipe for 
varying d w . This indicates that increasing roughness makes adsorption more difficult. 

Our results suggest an interesting scenario for the adsorption transition on fractal 
walls in more realistic models. First of all, for a continuous adsorption transition, 
4> is an increasing function of the wall dimension d w . Moreover, for high enough 
d w , (f> eventually reaches the value 1, its upper bound marking the onset of first-order 
adsorption. This fact is in open contrast with the results of Bouchaud and Vannimenus 
which, on the basis of scaling arguments, suggested the bounds 

l--^—(d L -d w )<cf><^. (14) 

u*saw 

We find that </> does not satisfy (|lj). Also the lower bound in ( |l4| ) is manifestly 
violated for low enough d w . Of course, here we deal with a hierarchical lattice, which 
is not fully adequate to represent consistently all the features of fractal objects. On 
the other hand, the bounds in ( |l4| ) were obtained by relying on formal analogies with 
cases of regular geometry. Similar violations of these bounds were reported previously, 
for both flat and fractal boundaries |)|, [l3| . Here we identify in the monotonicity of cf> 
and in the tendency of the transition to turn first-order as d w increases the physical 
reasons for the upper bound violation. In the investigation of Ref. E3 first-order 



7 



adsorption was found for a particular choice of lattice and fractal boundary among 
the many considered. A posteriori, we can understand that result as due to the fact 
that such choice determines a rather high d w relative to d saw . 

We obtained also results for random fractal walls. Unlike the case of deterministic 
walls above, this problem cannot be solved exactly, even on hierarchical lattices. 
However, quantities like (M) n and (N) n , which now have to be averaged also over 
wall randomness (overbar represents this average) , can be calculated quite accurately 
also for very large system sizes by a Monte Carlo approach 1 1 , 1^, |l8| . We performed 
these calculations for Cc with walls obtained by random combined applications of the 
rules in Figure 3c and 3d. At any level of lattice construction, we choose whether the 
wall is realized by rule 3c or 3d with probabilities 1 — A and A (A < 1), respectively. 
Of course, this determines which of the two generating functions, X 2 (Equation (12)) 
and A3 (Equation (13)), has to be used in order to calculate X n+ i in terms of X n . 
For < A < 1 the wall has, on average, a fractal dimension d w = ln(6 + 2A)/ln4. 
X n becomes now a random variable and we must consider its probability distribution. 
Of course, we can only produce a finite sampling of this distribution, by proceeding 
as follows. From a large set (up to 4 ■ 10 5 elements) {X n } of n-th level partition 
values we generate each element of the new sample {A n+ i} by choosing first, with the 
appropriate probabilities, between rules (12) and (13); then from {X n } are extracted 
at random the elements needed as entries into (12) or (13), and an element X n+ i of the 
new sample is computed as a function of them and of B n . Using some numerical tricks 
to control the possible rapid divergence of the partition functions near the transition, 
we could iterate this procedure up to n = 30 4- 35. 

By analyzing the scaling of (M) and (iV) n at the transition point (which now has 
to be numerically determined) as a function of L n , we could estimate 4> f° r different 
A's as reported in Figure 4. Even if the relative poorness of the samplings causes 
appreciable uncertainty in the <f> determinations (Figure 4), we see that this exponent 
stabilizes for A > 1 /2 to a value compatible with the upper limit <p — 1. The results are 
consistent with the scenario for deterministic walls. In particular, A = 1/2 corresponds 
to d w — 1.404, extremely close to d^ ~ 1.4 applying in that case. All this suggests 
that first-order adsorption should be expected also with random fractal walls, with 
the same threshold d* as in the deterministic case. 



3. Conclusions 



Our most remarkable result here is the roughness induced change into first-order of the 
adsorption transition. Without any other changes in the polymer-wall interactions, 
high enough d w make the adsorption transition discontinuous. This is found in all 
lattices considered, although marginally in Ca- The change in the nature of the 
transition takes place for d w definitely larger than d saw . In spite of the qualitative 
value of our model calculations, we can hope that similar properties could hold in 
realistic situations, in both 2D and 3D. 

The change into first-order of the adsorption transition should be imputed to the 
fact that, upon increasing d Wl the drop in entropy associated to a localization of the 
polymer near the wall, increases with d w . Even if it is not easy to give a precise 
meaning to the notion of distance of a polymer from a fractal substrate, at qualitative 
level we can argue that the entropic repulsion effect due to localization jl9|], should 
create a free energy barrier, whose height and (long) range certainly increase as d w gets 
larger. The rougher the surface, the more it limits the configurations of the confined 
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polymer. Thus, we can think d* w as the maximum dimension for which "tunneling" of 
the polymer can still occur continuously from the attractive free energy well at small 
distance to the unbound state at infinite distance across the barrier. For d w > d^, the 
"tunneling" becomes discontinuous because of the too large barrier (this means that 
the SAW, right at the desorption point, is not delocalized yet, as occurs in continuum 
adsorption) . 

A mechanism like the qualitative one outlined above has been demonstrated and 
precisely described for the phenomenon of wetting of self-affme rough substrates in 
2D, which has some analogies with our adsorption ]l2| |. Indeed, in that case one finds 
that a fluctuating interface depinns discontinuously from a rough substrate as soon 
as the roughness of the latter (measured by its self-affinity exponent £ w ) exceed £ , 
the exponent specifying the intrinsic roughness of the interface in the bulk (ll], |l2|| . 
The coincidence of the threshold roughness with (o (which by analogy, would suggest 
dw = here) is a peculiar feature of the interfacial problem in 2D. Indeed, for that 
problem, the possibility of a path-integral description in 1 + 1 dimensions allows to 
establish a correspondence with the ID quantum tunneling of a particle across a long 
range repulsive potential barrier. It turns then out that a roughness £ w — (o = 1/2 
determines a decay of this potential right at the threshold for discontinuous tunneling 
p2{ . Our results here show that fractal wall roughness leaves room for a continuous 
polymer adsorption also when d w > d saw . 
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Table 1. SAW adsorption on deterministic fractal walls in C\. Wall rules are applied 
in the reported order. 



d L = In 5/ In 2 = 2.3219, 


u> c = 0.3660, d saw 
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X 2 , X\ 
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x 2 


1.5850 
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Table 2. SAW adsorption on deterministic fractal walls in Cb- 



d L = 3, u> c = 


0.2441, d saw = 


1.1995 
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X! 
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3.3049 
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5/4 
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X2 , X\ , X\ 


4/3 


3.3500 
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3/2 


3.3549 


0.9637 
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3.3629 


l a 


X2.X2.X2, X\ 
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3.3652 
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3.3663 
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a Ai > A s 



Table 3. SAW adsorption on deterministic fractal walls in Cq- 



d L = In 12/ 


ln4 = I.' 


7925, u) c = 


0.5175, dsa-w = 1.0649 




d w 


kc 


4> 


Xi 


1 


1.3279 


0.5437 


x 2 


1.2925 


1.4715 


0.8140 


X3,X2 


1.3962 
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a Ai > A s 



10 




Figure 1. a Construction rule of Ca- The length of the diagonal step is different 
from the other lattice steps only for drawing convenience, b-d Construction rules of 
deterministic fractal walls (heavy lines) in Ca- The allowed polymer configurations 
develop to the right of the walls. 




Figure 2. a Construction rule of Cb- b—c Construction rules of deterministic fractal 
walls (heavy lines) in Cb- 




Figure 4. <j> as a function of A for SAW adsorption on random fractal walls in Cc- 
4> values at A = and A = 1 (closed circles) are exact. 



